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Abstract 

Elliptical galaxies are modelled as homeoidally striated Jacobi ellip- 
soids (Caimmi & Marmo 2005) where the peculiar velocity distribu- 
tion is anisotropic, or equivalently as their adjoints configurations i.e. 
classical Jacobi ellipsoids of equal mass and axes, in real or imag- 
inary rotation (Caimmi 2006). Reasons for the coincidence of bi- 
furcation points from axisymmetric to triaxial configurations in both 
the sequences (Caimmi 2006), contrary to earlier findings (Wiegandt 
1982a, b; Caimmi & Marmo 2005) are presented and discussed. The 
effect of centrifugal support at the ends of the major equatorial axis, 
is briefly outlined. The existence of a lower limit to the flattening of 
elliptical galaxies is investigated in dealing with a number of limiting 
situations. More specifically, (i) elliptical galaxies are considered as 
isolated systems, and an allowed region within Ellipsoidland (Hunter 
& de Zeeuw 1997), related to the occurrence of bifurcation points from 
ellipsoidal to pear-shaped configurations, is shown to be consistent 
with observations; (ii) elliptical galaxies are considered as embedded 
within dark matter haloes and, under reasonable assumptions, it is 
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shown that tidal effects from hosting haloes have little influence on 
the above mentioned results; (iii) dark matter haloes and embedded 
elliptical galaxies, idealized as a single homeoidally striated Jacobi el- 
lipsoid, are considered in connection with the cosmological transition 
from expansion to relaxation, by generalizing an earlier model (Thuan 
& Gott 1975), and the existence of a lower limit to the flattening of 
relaxed (oblate-like) configurations, is established. On the other hand, 
no lower limit is found to the elongation of relaxed (prolate-like) con- 
figurations, and the observed lack of elliptical galaxies more elongated 
than E7 needs a different physical interpretation, such as the occur- 
rence of bending instabilities (Polyachenko k, Shukhman 1979; Merritt 
& Hernquist 1991). 

keywords - cosmology: dark matter - galaxies: evolution - galaxies: 
formation - galaxies: haloes - galaxies: structure. 

1 Introduction 

Large-scale celestial objects, represented as self-gravitating fluids, exhibit 
different features according if their subunits, or "particles" , are conceived as 
"collisional" or "collisionless" . 

In the former alternative, the gravitational field that is generated by the 
system as a whole is negligible in respect of the force between two colliding 
subunits, when they repel each other strongly. Consequently, particles in 
self-gravitating, collisional fluids, are subjected to violent and short-lived 
accelerations as they are sufficiently close to each other, interspersed with 
longer periods when they move at nearly constant velocity. 

In the latter alternative, the gravitational field that is generated by the 
system as a whole is dominant in respect of the force between two sub- 
units, even if they are close to each other. Consequently, particles in self- 
gravitating, collisionless fluids, are subjected to smooth and long-lived accel- 
erations through their trajectories. 

Gas in stars and stars in stellar systems may be approximated, to a good 
extent, as collisional and collisionless, ideal self-gravitating fluids, respec- 
tively. The statistical problem of particle motion in the latter case is that of 
the former, with the collisions left out. "Ideal" has to be intended as "par- 
ticles collide as perfectly and undeformable spheres" in the former situation, 
and "particles do not interact each other at all" in the latter. 
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The motion equation of fluid flow turns out to be the same for both 
collisional and collisionless fluids (e.g., Jeans 1929, Chap. II, §26, Chap. VII, 
§§ 211-215; Binney & Tremaine 1987, Chap. 4, § 4.2) provided gases are thought 
of as far from equilibrium. Then the velocity distribution of molecules no 
longer obeys the Maxwell distribution and is, in general, anisotropic; accord- 
ingly, the pressure is represented by a stress tensor. In the special case of 
isotropic velocity distributions, the pressure attains its usual meaning and 
the motion equation reduces to the Euler's equation, provided the velocity 
of an infinitesimal fluid element is intended as the mean velocity of all the 
particles within the same element at the time considered. On the contrary, 
the Euler's equation for a given fluid with isotropic distribution can be gen- 
eralized to an anisotropic one, by replacing velocities with mean velocities 
and pressures with stress tensors, in the sense specified above. 

For special classes of ideal, self-gravitating fluids, such as steadily ro- 
tating polytropes with polytropic index n > 1/2 (Vandervoort 1980a) and 
isothermal spheres (e.g., Binney & Tremaine 1987, Chap. 4, §4.4b), a one-to- 
one correspondence has been discovered between collisional and collisionless 
systems with equal physical parameters. Steadily rotating polytropes may 
be represented, to a first extent, as steadily rotating, homeoidally striated 
ellipsoids (e.g., Vandervoort 1980b; Vandervoort & Welty 1981; Lai et al. 
1993). 

Though most astronomical bodies exhibit ellipsoidal-like shapes and isopy- 
cnic surfaces, still some caution must be used in dealing with the above 
approximation, with regard to local values of physical parameters. This 
is why steadily rotating polytropes are in hydrostatic equilibrium while, 
owing to the Hamy's theorem (e.g., Chambat 1994), the contrary holds 
for their homeoidally striated ellipsoidal counterparts. On the other hand, 
homeoidally striated ellipsoids may safely be assumed as a viable approx- 
imation to self-gravitating fluids, with regard to typical values of physical 
parameters, averaged over the whole volume 1 (e.g., Vandervoort 1980b; Van- 
dervoort & Welty 1981; Lai et al. 1993). 

Large-scale, celestial objects have been modelled as collisional, self-gravitating 
fluids, in particular homeoidally striated ellipsoids, since the beginning of 
their classification (see e.g., Chandrasekhar 1969, Chap. 1). The evidence 

1 An important exception is the energy ratio of rotational to random motions, E rot /E pec 
(Caimmi 1979, 1983). 
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of rotation in spiral galaxies and the symmetry of figure shown by ellip- 
tical galaxies and spiral bulges, suggested the following (e.g., Jeans 1929, 
Chap. XIII, §299): (i) all (regular) galaxies rotate; (ii) the symmetry of 
figure is precisely such as rotation might be expected to produce; (iii) the 
observed shapes of (regular) galaxies can be explained as the figures assumed 
by masses rotating under their own gravitation. In particular, the system of 
primeval stars must have conserved roughly the same shape as the original 
mass of gas from which it emerged (e.g., Blaauw 1965). The lack of elliptical 
galaxies more flattened than E7 was explained in a classical paper (Thuan & 
Gott 1975) where virialized elliptical galaxies and their parent density per- 
turbations are modelled as MacLaurin spheroids and homogeneous, rigidly 
rotating spheres, respectively. 

Since then, observations began to yield increasing evidence that (giant) 
elliptical galaxies cannot be sustained by systematic rotation (e.g., Bertola 
& Capaccioli 1975; Binney 1976; Illingworth 1977, 1981; Schechter & Gunn 
1979). Accordingly, (giant) elliptical galaxies were conceived as collisionless, 
self-gravitating fluids, with triaxial 2 boundaries set by specific anisotropic pe- 
culiar velocity distribution of stars (Binney 1976, 1978, 1980), and a negligi- 
ble contribution from figure rotation. Owing to high-resolution simulations, 
the same holds also for (nonbaryonic) dark matter haloes hosting galaxies 
and cluster of galaxies (e.g., Hoeft et al. 2004; Rasia et al. 2004; Bailin & 
Steinmetz 2004). 

Isotropic peculiar velocity distributions necessarily imply configurations 
which rotate around the minor axis. Accordingly, empirical evidence of sys- 
tematic rotation around the major axis (e.g., Bertola & Galletta 1978), in 
absence of tidal potential, makes a signature of the occurrence of anisotropic 
peculiar velocity distribution. 

Though anisotropy in peculiar velocity distribution is a basic ingredient in 
the description and investigation of stellar systems and hosting dark matter 
haloes, still no attempt (to the knowledge of the author) has been made to 
explain why elliptical galaxies more flattened than E7 do not exist, following 
the same line of thought as in Thuan & Gott (1975), with regard to collision- 
less fluids. This paper aims to address this lack, and is based on a theory 

2 Strictly speaking, all tridimensional bodies may be conceived as triaxial, in particular 
spheres and spheroids. Throughout this paper, "triaxial" shall be intended as denoting 
ellipsoids where the axes are different in length. 
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that systematic and random motions are unified (Caimmi 2006, hereafter 
quoted as C06 3 ). To this aim, the procedure used by Thuan & Gott (1975) 
shall be extended from MacLaurin spheroids to homeoidally striated Jacobi 
ellipsoids (Caimmi & Marmo 2005, hereafter quoted as CM05 4 ), where the 
effect of peculiar motion excess, implying anisotropic peculiar velocity distri- 
bution, is equivalent to the effect of additional (real or imaginary) rotation, 
related to isotropic peculiar velocity distribution (C06). 

An alternative explanation for the absence of elliptical galaxies (and non- 
baryonic dark matter haloes) more flattened or elongated than E7 as due 
to bending instabilities, has been suggested long time ago from analytical 
considerations involving homogeneous (oblate and prolate) spheroids (Poly- 
achenko & Shukhman 1979; Fridman & Polyachenko 1984, Vol.1, Chap. 4, 
Sect. 3.3, see also pp. 313-322; Vol. 2, p. 159) and numerical simulations in- 
volving inhomogeneous (oblate and prolate) spheroids (Merritt & Hernquist 
1991; Merritt & Sellwood 1994). 

In a cosmological scenario (Thuan & Gott 1975), the occurrence of a 
limiting ellipticity in oblate configurations depends on the amount of spin 
growth regardless from the onset of bending instabilities. An interesting 
question could be if a similar conclusion holds for prolate configurations. 
It will be seen that the onset of bending instabilities is necessary for the 
occurrence of a limiting ellipticity in prolate configurations. 

This paper is structured in the following manner. The general theory of 
homeoidally striated Jacobi ellipsoids, and the relevant results, are outlined 
in Sect. 2. A unified theory of systematic and random motions, implying 
a definition of imaginary rotation, and the relevant results, are outlined in 
Sect. 3. An extension of Thuan & Gott (1975) procedure to homeoidally stri- 
ated Jacobi ellipsoids in real or imaginary rotation, and the relevant results, 
are outlined in Sect. 4. The existence of a lower limit to the flattening of 
elliptical galaxies is investigated in Sect. 5 taking into consideration a num- 
ber of simplified situations, namely (i) elliptical galaxies as isolated systems; 
(ii) elliptical galaxies as embedded in dark matter haloes; (iii) dark matter 
haloes and hosted elliptical galaxies, idealized as a single homeoidally stri- 
ated Jacobi ellipsoid, in connection with the cosmological transition from 

3 A more extended file including an earlier version of the current paper is available at 
the arxiv electronic site, as astro-ph/0507314. 

4 A more extended version is available at the arxiv electronic site, as astro-ph/0505306. 
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expansion to relaxation. Some concluding remarks are drawn in Sect. 6, and 
a few arguments are treated with more detail in the Appendix. 



2 General theory 

2.1 Homeoidally striated Jacobi ellipsoids 

A general theory for homeoidally striated density profiles has been developed 
in earlier approaches (Roberts 1962; Caimmi 1993a; Caimmi & Marmo 2003, 
hereafter quoted as CM03; CM05), and an interested reader is addressed 
therein for deeper insight. What is relevant for the current investigation, 
shall be mentioned and further developed here. 

The isopycnic (i.e. constant density) surfaces are defined by the following 
law: 



where p , r , are a scaling density and a scaling radius, respectively, related 
to a reference isopycnic surface, £ is a scaled distance, independent of the 
direction along which radial coordinates are calculated, and S = £(r), is 
related to the boundary, where r = R. Both cored and cuspy density profiles, 
according to the explicit expression chosen for the scaled density, /(£), are 
represented by Eqs. (1). For the cored profiles, a different normalization is 
used here with respect to Caimmi (1993a), where £ = r/R and p is the 
central density. 

The mass, the inertia tensor, and the potential self-energy tensor are: 



P = Pof(0 




/(1) = 1 ; Po = p(1) 



(la) 
(lb) 



M = z/ mas M ; 

Ipq = S pq U inr Mal 



(2) 
(3) 




(4) 



(6) 



(5) 
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<S pq — u se i(B se i) pg ; S — u se iB se \ ; (7) 

where 5 pq is the Kronecker symbol; G is the constant of gravitation; z/ mas , 
^mr, ^sei, are profile factors i.e. depend only on the density profile via the 
scaled radius, S; a±, a 2 , a 3 , are semiaxes; e pq = a p /a q are axis ratios; Ai, A2, 
A 3 , are shape factors i.e. depend only on the axis ratios; and M is the mass 
of a homogeneous ellipsoid with same density and boundary as the reference 
isopycnic surface: 

An 

M = —poaoia 02 a 03 ; (8) 

where a i, a 02 , a 03 , are the semiaxes of the ellipsoid bounded by the reference 
isopycnic surface. The combination of Eqs. (lb), (2), and (8) yields: 



P ^mas 

where p = 3M/ (Ananas) is the mean density of the ellipsoid. 
The limiting case of homogeneous configurations reads: 



(9) 



/(0 = 1 , 0<£<H ; (10a) 

1 3 

z/ mas = H 3 ; z/ inr = - ; u sel = — ; (10b) 
10 

for further details, see CM03, CM05. 
2.2 Systematic rotation 

In dealing with angular momentum and rotational energy, the preservation 
of (triaxial) ellipsoidal shape imposes severe constraints on the rotational 
velocity field. In the special case of homeoidally striated Jacobi ellipsoids, 
the systematic velocity field is defined by the law (CM05): 

v Iot (r,9,<f>) = ^(01,71-/2,0) 

v rot (R,6,(j)) ~ Urot (ai,7r/2,0) ' 1 j 

or equivalently, the angular (with respect to x 3 axis) velocity field is defined 
by the law: 

^K,7r/2,0) 

n(R,0,(f)) fi(ai,7r/2,0) ' 1 ' 
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where (r, 9, <p), (R, 6, 0), represent a point on a generic isopycnic surface and 
on the boundary, respectively, along a fixed radial direction, and (a[, ir/2, 0), 
(ai,7r/2,0), represent the end of the major equatorial semiaxis. 

It is worth noticing that the following special cases are described by 
Eqs. (11) or (12): (a) rigid rotation; (b) constant rotation velocity every- 
where; (c) rigid rotation of isopycnic surfaces and constant rotation velocity 
on the equatorial plane. To maintain ellipsoidal shapes, differential rota- 
tion [e.g., cases (b) and (c) above] must necessarily be restricted to axysim- 
metric figures i.e. spheroidal configurations. In the limiting situation of 
homogeneous, rigidly rotating, dynamical (or hydrostatic) equilibrium con- 
figurations, homeoidally striated Jacobi ellipsoids reduce to classical Jacobi 
ellipsoids (CM05). 

The angular-momentum vector and the rotational-energy tensor are: 

Js = 5 S 3?7anm^anmMap(l + e 2 qp )(v lot ) p ] p ± q ± S ] (13) 

{E rot )p q = 5 pq (l - 5 p3 )r] ro tU TOt M {(v TOt ) p } 2 ; (14) 
and the related module and trace, respectively, read: 

J = ??anm^anmM(l + C^Ol (f rot) 1 5 (15) 
^rot = ^rot^rotM(l + 6^) [(w rot )l] 2 5 (16) 

where r) anm , Vrot, are shape factors, u amni z/ rot , are profile factors, (v TOt ) p is the 
rotational velocity at the end of the semiaxis, a p , p = 1,2, and the rotation 
axis has been chosen to be x 3 . For further details, see CM03, CM05. 
The combination of Eqs. (15) and (16) yields: 

J 2 J 2 

(E rot ) pq = ^ram {B Taxn ) p g = j^^2 ^P<? ' (17a) 

^ram = \ (17b) 



Z/ 2 

anm 



•2 



(5 ram )p 9 = 5 m {\ - 5 p3 )^^ pl ; (17c) 

Vanm I 1 + e 2l) 

T^"pq ^ram(-Sram)p9 i (17d) 

which makes an alternative expression of the rotation-energy tensor, and: 

J 2 J 2 

-E'rot — o^ramSram = — 2 ', (18a) 

Ma{ Maf 
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^anm 1 + e 21 

^ = ^ram-Bram ! (18°) 

which makes an alternative expression of the rotation energy. 
The limiting situations outlined above, read: 

^anm = ^rot = ^inr = ~ \ ^ram = ^inr = 5 ; (19) 



for homogeneous configurations in rigid rotation, according to case (a); 

1 1 16 

^anm — ^ j ^rot — g j ^ram — j v^UJ 

for homogeneous configurations with constant velocity on the equatorial 
plane, according to cases (b) and (c); 

1 . . 

?7anm = 1 5 ^rot = ~ 5 (21) 

for rigidly rotating isopycnic surfaces, cases (a) and (c); 

37T 3 , . 

for constant velocity everywhere, case (b). Further details can be found in 
CM03, CM05. 

2.3 Virial equilibrium configurations 

Let us define virial equilibrium as characterized by the validity of the virial 
theorem, and relaxed and unrelaxed configurations as systems where virial 
equilibrium does and does not coincide, respectively, with dynamical (or 
hydrostatic) equilibrium (CM05). 

With regard to unrelaxed configurations, the generalized tensor virial 
equations read (CM05): 

(E S e\) P q + 2(-E ro t) w + 2( pq E pec = ; (23) 
C P9 = %^ ; P= 1,2,3 ; q = 1,2,3 ; (24) 

EC PP = f^ = C ; o<c PP <C ; (25) 

p=l ^pec 

C P9 = ; p^q ; (26) 
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where ( pq is the generalized anisotropy tensor, E pcc is the residual energy, 
and (E pcc ) pq = ( pq E pcc is the effective residual-energy tensor i.e. the right 
amount needed for the configuration of interest to be relaxed. The diagonal 
components of the generalized anisotropy tensor, ( pp , may be conceived as 
generalized anisotropy parameters. The related trace, (, may be conceived 
as a virial index, where ( = 1 corresponds to null virial excess, 2AE pec = 
2(E pcc — -Ep CC ), which does not necessarily imply a relaxed configuration 5 , 
C > 1 to positive virial excess, and ( < 1 to negative virial excess. 

The substitution of Eqs. (4)-(7) and (17)-(18) into the virial equations 
(23), and the particularization to the rotation axis, p = 3, allows the following 
expression of the peculiar energy (CM05): 

E pcc - - — — — , {!() 
1 «1 (,33 

accordingly, the remaining virial equations read (CM05): 

{( 3 3S qq - ( qq S 33 ) - 2h( 33 n qq = ; g = l,2 ; (28) 

where h may be conceived as a rotation parameter. It is apparent that 
Eq. (28) admits real solutions provided the following inequalities hold: 

k<|i . g = l,2 ; (30) 

S33 <-"33 

which is the natural extension of its counterpart related to axisymmetric, 
relaxed configurations (Wiegandtl982a,b). 

2.4 Rotation and anisotropy parameters 

To get a more evident connection with the centrifugal potential on the bound- 
ary, let us define the rotation parameter (CM05): 

" = ; (31) 



5 For instance, a homogeneous sphere undergoing coherent oscillations exhibits £ > 1 
at expansion turnover and £ < 1 at contraction turnover. Then it necessarily exists a 
configuration where £ = 1 which, on the other hand, is unrelaxed. 
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where Q = Q(ai, 0, 0) is the angular velocity at the end of the major equato- 
rial semiaxis, denoted as ai, and p is the mean density of the ellipsoid: 

3 M 

p = - A ; 32 

the above definition of the rotation parameter, v, makes a generalization of 
some special cases mentioned in the literature (e.g., Jeans 1929, Chap. IX, 
§232; Chandrasekhar & Leboviz 1962). 

The combination of Eqs. (15), (29), (31), and (32) yields: 

i 3 2 2 (1 e 2l) 2 / 00 \ 
h = ^anm^m— V ] (33) 

which links the rotation parameters, h and v. An explicit expression of the 
rotation parameter, h, may directly be obtained from Eq. (28), as: 

h= 1 ^-^ ; , = i,2 ; (34) 

* S33 Kqq 

and the substitution of Eqs. (34) into (33), using (17) and (18), allows the 
explicit expression of the rotation parameter, v, as: 



^U + ^lJ Vamn^amn S33 '<-gg 



which, in the special case of rigidly rotating, homogeneous configurations 
with isotropic peculiar velocity distribution, reduces to a known relation for 
Jacobi ellipsoids and, with the additional demand of axial symmetry, to a 
known relation for MacLaurin spheroids (e.g., Jeans 1929, Chap. VIII, §§189- 
193; Chandrasekhar 1969, Chap. 5, §32, Chap. 6, §39; Caimmi 1996a). 

An explicit expression of anisotropy parameter ratio, can be obtained via 
Eq. (34). The result is (CM05): 



1Z 

S qq 



Cqq_ _ Sqq 
C33 £33 

Cn _ gn - 2hn n 

C22 <^22 — 2/i7?-22 



9 = 1,2 ; (36) 

(37) 
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and the combination of Eqs. (25) and (36) yields: 



C33 £33 



(38) 



C S-2KR ' 

which provides an alternative expression to Eqs. (34) and (35), as (C06): 

fr = - C33tS ~ CtS33 ; (39) 

_ 1 £21 £31 1 C33i> ~ C^33 

v ~ Hi + e 2 21 y V i mn ^ amn ( 33 n ; 1 } 

that are equivalent to Eq. (28), and then admit real solutions provided in- 
equality (30) is satisfied. 

Finally, Eqs. (28) may be combined as: 

T^ll C33<Sll — Cll^33 . 



C33^22 — C22<5; 



33 



where it can be seen that Eqs. (37) and (41) are changed one into the other, 
by replacing the terms, SzzCqq/C^i with the terms, 2h1Z qq , and vice versa. 

In the special case of axisymmetric configurations, the shape factors re- 
lated to equatorial axes coincide, A 2 — A 1: as e 2 i = 1 (e.g., CM05). Then 
<Sn = S22 owing to Eqs. (6), (7), and Tin = IZ22 owing to Eqs. (17), which 
necessarily imply (u = (22, owing to Eq. (41). 

In the general case of triaxial configurations, the contrary holds, A2 ^ Ai, 
as €21 7^ 1, then Su 7^ S22 and IZu 7^ ^22, then the equality, C11 = C22, via 
Eq. (37), implies the validity of the relation: 

1 S 11 -S 22 

2nn-n 2 2 1 ; 

if otherwise, the residual velocity distribution along the equatorial plane 6 is 
anisotropic i.e. C11 (22- The related degeneracy can be removed using an 
additional condition, as it will be shown in the next section. 



6 Throughout this paper, "along the equatorial plane" has to be intended as "along any 
direction parallel to the equatorial plane" . 
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3 A unified theory of systematic and random 
motions 



3.1 Imaginary rotation 

A unified theory of systematic and random motions is allowed, taking into 
consideration imaginary rotation, as discussed in earlier attempts (Caimmi 
1993b; C06), and an interested reader is addressed therein for deeper insight. 
What is relevant for the current investigation, shall be mentioned and further 
developed here. 

It has been shown above that Eq. (28), or equivalently one among (34), 
(35), (39), (40), admits real solutions provided inequality (30) is satisfied. If 
otherwise, the rotation parameter - let it be h or v - has necessarily to be 
negative, which implies, via Eqs. (29) or (31), an imaginary angular velocity, 
ifl, where i is the imaginary unit. Accordingly, the centrifugal potential takes 
the general expression (C06): 

T T (xi,x 2 ,x 3 ) = T^[^(x 1 ,x 2 ,x 3 )] 2 w 2 ; w 2 = x\ + x\ ; (43) 

where the minus and the plus correspond to imaginary and real rotation, 
respectively. The centrifugal force related to real rotation, dT + /dw, has 
opposite sign with respect to the gravitational force, dV/dw. On the other 
hand, the centrifugal force related to imaginary rotation, dT~ /dw, has equal 
sign with respect to the gravitational force, dV/dw. Then the net effect 
of real rotation is flattening, while the net effect of imaginary rotation is 
elongation, with respect to the rotation axis (Caimmi 1996b). 

To get further insight, let us particularize Eq. (28) to the special case of 
null rotation (h — 0). The result is: 

= 5 h = ; 9 = 1,2 ; (44) 

S33 «->33 

where the right-hand side, via Eqs. (6) and (7), depends on the axis ratios 
only. It can be seen (Appendix A) that S qq /S 3 3 > 1 for oblate-like con- 
figurations (ai > a 2 > £13), and S qq /S 33 < 1 for prolate-like configurations 
(a 2 < ai < a 3 ), which implies C99/C33 ^ 1 f° r oblate-like configurations, 
and C99/C33 < 1 f° r prolate-like configurations. Accordingly, the net effect 
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of positive or negative residual motion excess along the equatorial plane is 
flattening or elongation, respectively. In what follows, it shall be intended 
that residual motion excess is related to the equatorial plane. 



3.2 Residual motion excess and rotation 



With regard to the rotation parameter, v , it is convenient to use the more 
compact notation 7 (C06): 



S^rot^rot 

v N = V 

^sel 



(45) 



and the substitution of Eqs. (6), (7), (17), (18), into (35) and (40), yields the 
equivalent expressions: 



1 



9=1,2 ; 



1 + e 2 



21 



2 i , C33 



A 1 + etA 2 + 



C: 



4 

— e 3l A 3 



33 



(46) 
(47) 



which shows that, for homeoidally striated Jacobi ellipsoids, the normalized 
rotation parameter, t>N, depends on the axis ratios, 621, £31, and a single 
anisotropy parameter, £33 = C33/C 111 the special case of homogeneous, 
rigidly rotating configurations, owing to Eqs. (10b), (19), and (21), Eq. (45) 
reduces to: = v. 

In the limit of isotropic residual velocity distribution, (u = (22 = C33 = 
C/3, Eqs. (46) and (47) take the simpler form (C06): 



Mi 
Mi 



1,2 



1 



1 + & 



A 1 + e 2 21 A 2 - 2e 2 31 A 3 



(48) 
(49) 



where the index, iso, means isotropic residual velocity distribution. 
Accordingly, Eqs. (46) and (47) may be expressed as (C06): 

= (^N)iso - (^N)ani 5 



(50) 



7 The factor, 3, does not appear in CM05: wn = 3(un)cmo5- The current normalization 
is more convenient in dealing with a general theory (Caimmi, in preparation) . 
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Ma„i = h« - 1 elA 3 ; (51) 




Mani = ( 7^ - 3 ] ; (52) 

where (i>N)am > for oblate-like configurations, Cqq/Cs3 > lj ( w N)am < for 
prolate-like configurations, ( qq /(33 < 1; and the index, ani, means contri- 
bution from residual motion excess. Accordingly, positive residual motion 
excess is related to real rotation. On the contrary, negative residual motion 
excess is related to imaginary rotation. 
Let us rewrite Eq. (50) as: 

(^N)iso = + (^N)ani 5 (53) 

which, owing to Eqs. (31) and (45), is equivalent to: 

nl = n 2 T ^L ; (54) 

where the plus corresponds to real rotation, (t>N)am > 0, and the minus to 
imaginary rotation, (fN)ani < 0. Then the effect of residual motion excess 
on the shape of the system, is virtually indistinguishible from the effect of 
additional systematic rotation, with similar (rotation) velocity distribution as 
the pre-existing one. The related, explicit expression, according to Eqs. (12), 
(31), and (52), is (C06): 

flaniM) = ^aniK,0) = »K,0) = gM) ( 

«L = nL(ai) = T{v^U^^2nG-p ; (56) 

J^rot^rot 

where the double sign ensures a non negative value on the right hand-side 
member of Eq. (56). Accordingly, a homeoidally striated Jacobi ellipsoid with 
assigned systematic rotation and residual velocity distribution, with regard 
to the shape, is virtually indistinguishable from an adjoint configuration of 
equal mass and axes, systematic rotation velocity distribution deduced from 
Eqs. (12), (55) and (56), and isotropic residual velocity distribution. 

Owing to Eqs. (50), (51), and (52), Eqs. (48) and (49) are equivalent to 
(46) and (47), respectively. On the other hand, Eqs. (48) and (49) are valid 
for classical Jacobi ellipsoids with rotation parameter, v = (v^)i SO - 
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Accordingly, a homeoidally striated Jacobi ellipsoid (in general, an unre- 
laxed configuration) with assigned systematic rotation and residual velocity 
distribution, with regard to the shape, is virtually indistinguishable from an 
adjoint, classical Jacobi ellipsoid (a relaxed configuration) of equal mass and 
axes, and rotation parameter equal to the normalized rotation parameter of 
the original configuration. 



3.3 Axis ratios and anisotropy parameters 

The combination of the alternative expressions of the rotation parameter, 
(^N)iso, expressed by Eqs. (48), yields: 

4(A 2 - Ax) = (l-4)4As ; (57) 

which, for axisymmetric configurations, reduces to an indeterminate form, 
= 0. 

The combination of the alternative expressions of the rotation parameter, 
(^N)ani, expressed by Eqs. (51), yields: 

C33-C22 = e!i(C33-Cll) ; (58) 

which, for isotropic residual velocity distributions, reduces to an indetermi- 
nate form, = 0. In addition, axisymmetric configurations (e 2 i = 1) nec- 
essarily imply isotropic residual velocity distributions along the equatorial 
plane, Cn = (22- 

The combination of Eqs. (25) and (58) yields: 



C - (2 - 4)0 



C11 = - r — ^ 33 ; (59a) 

1 T e 21 



4C + (1-24)0 



C22 = — — 2 , (59b) 

1 T e 21 

which, for axisymmetric configurations (e 2 i = 1) reduces to £n = £22 = 
(C - C33V2, and the special case, C 33 = C/3, reads C11 = C22 = C/3- 

The condition, (u > 0, related to Eqs. (24), necessarily implies < 
C33/C ^ 1/2 for e 2 i = 0. If otherwise, sequences of virial equilibrium configu- 
rations cannot attain the oblong shape, e 2 i = e 3i = 0, but must stop earlier, 
where the equivalent relations: 

e 2 i = — 7 I C11 = ; (60) 

<533 
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are satisfied. 

Accordingly, with regard to homeoidally striated Jacobi ellipsoids, the 
anisotropy parameters along the equatorial plane, Cn an d C22) cannot be 
arbitrarily assigned, but depend on the equatorial axis ratio, e 2 i, conform 
to Eqs. (59). On the other hand, the further knowledge of the meridional 
axis ratio, 631, and the rotation parameter, v^, allows the determination of 
the rotation parameter, (i^ani, via Eqs. (48), (49), (50), and then the ratios, 
Cqq/C^) C33/C) y i a Eqs. (51), (52), respectively, or the anisotropy parameter 
along the rotation axis, C33, provided the virial index, (, defined by Eq. (25), 
is assigned. 

In conclusion, with regard to homeoidally striated Jacobi ellipsoids de- 
fined by assigned axis ratios, €31, €21, rotation parameter, v^, and virial index, 
(, the anisotropy parameters, Cii> (22, (33, cannot arbitrarily be fixed, but 
must be determined as shown above. 

3.4 Centrifugal support along the major equatorial axis 

The calculation of the gravitational force, induced by homeoidally striated 
Jacobi ellipsoids, involves numerical integrations (e.g., Chandrasekhar 1969, 
Chap. 3, § 20) and is outside the aim of the current paper. With regard to an 
assigned isopycnic surface internal to, or coinciding with, the boundary, the 
gravitational force, F G , induced at the end of the related major equatorial 
semiaxis, a[, satisfies the inequality: 

-[F G (a[, 0, 0)] sph < -F G (a[, 0, 0) < -[F G (a[, 0, 0)] foc ; 

where the left-hand side is related to spherical isopycnics with unchanged 
major equatorial axes, and the right-hand side to confocal isopycnic surfaces 
from the one under consideration to the centre. Owing to the Newton's and 
MacLaurin's theorem (e.g., Caimmi 2003), the following relations hold: 

F G (a[,0,0) = -27rGp xxx (a / 1 )(A 1 ) xxx a , 1 ; (61) 

where p(a[) is the mean density within the isopycnic surface, and xxx = sph, foe, 
for the striated sphere and the focaloidally striated ellipsoid surrounded by 
the homeoidally striated corona, respectively. 

The balance between gravitational and centrifugal force at the end of the 
major equatorial axis of the isopycnic surface under discussion, reads: 

-27rGp xxx (a / 1 )(A 1 ) xxx a , 1 + J&KK = 5 ( 62 ) 
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where the angular velocity, ft- lso , takes into account both systematic rotation 
and residual motion excess, according to Eq. (54). On the other hand, the 
generalization of Eq. (31) to a generic isopycnic surface within the boundary, 
reads: 

^° = 9 J „ M ; (63) 

and the combination of Eqs. (45), (53), (54), (62), (63), yields: 

{Ko«)Uxxx = 9 = (^i)xxx ; (64) 

where the index, eq, denotes the rotation parameter, fi so (ai), related to 
centrifugal support at the end of major equatorial axis of the isopycnic surface 
under consideration. The above results allow the validity of the relation: 

{Ko(ai)]eq}cof < bisoK)]eq < { KoK )]cq}sph \ (65) 

where [fiso(oi)]eq is the critical value related to the homeoidally striated 
Jacobi ellipsoid, with regard to an assigned isopycnic surface internal to, 
or coinciding with, the boundary. Then v- lso {d^) > {[fi S o( a i)]eq}cof and 
^iso( a 'i) < {[^iso( a 'i)]cq}sph make a sufficient and a necessary condition, re- 
spectively, for the occurrence of centrifugal support at the end of major 
equatorial axis in the case under consideration. 

Owing to the Newton's and MacLaurin's theorem, the condition of cen- 
trifugal support, Eq. (64), for focaloidally striated ellipsoids surrounded by 
homeoidally striated coronae, coincides with its counterpart related to ho- 
mogeneous ellipsoids with equal mean density, axis ratios, and velocity field. 
On the other hand, the rotation parameter, fi S o(a'i), in the special case of 
focaloidally striated ellipsoids, is also expressed by Eq. (64), particularized to 
classical Jacobi ellipsoids. A comparison with the last part of Eq. (64) shows 
that centrifugal support in focaloidally striated ellipsoids occurs only for flat 
configurations, e 3 i = 0. Accordingly, Eq. (65) reduces to: 

2 

< KoK)]cq < - ; (66) 

owing to Eq. (64) and (Ai) sph = 2/3 (e.g., CM03). 

With regard to rigid rotation, it is a well known result that the occurrence 
of centrifugal support, at the end of the major equatorial axis, depends on 
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the steepness of the density profile (e.g., Jeans 1929, Chap. IX, §§230-240). 
More precisely, a steeper density profile implies an earlier centrifugal support 
and vice versa, unless the bifurcation point from ellipsoidal to pear-shaped 
configurations is attained, which occurs for nearly homogeneous matter dis- 
tributions. A similar trend is expected to hold for any velocity profile of the 
kind considered in the current paper. 

The combination of Eq. (65) with its counterpart related to the boundary, 
yields: 

M fa'A 3 ft 2 (a'i) 

a further restriction to mass distributions obeying the following law: 
M(a') (a'\ k 

where k = 3,1,0, represent homogeneous, isothermal, and Roche-like mass 
distributions, make Eq. (67) reduce to: 

vM =^o - 7^p{ ; (69) 

\aij i^(ai) 

where the last factor equals unity for rigid rotation and the ratio, (ai/a^) 2 , for 
constant velocity along the major equatorial axis. Accordingly, centrifugal 
support at the end of major equatorial axis is first attained on the boundary 
in the former alternative, fi so (a'i) < ^iso(oi)- 
In the latter alternative, Eq. (69) reduces to: 

( a ' Y~ k 

^iso(ai) = v iso ( -ij ; (70) 

which, for sufficiently steep density profiles, < k < 1, shows a similar trend 
with respect to the former alternative. On the other hand, centrifugal support 
is first attained everywhere along the major equatorial axis for isothermal 
mass distributions, k — 1, and at the centre for sufficiently mild density 
profiles, 1 < k < 3, which implies Vi SO (a[) > Vi SO (ai). 
The combination of Eqs. (45) and (64) yields: 

{[(^N)iso(a / l)]cq}xxx = ^^^ (^Oxxx 5 (71) 

^sel 
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and Eq. (66) reads: 



o < [MisoK)] eq < ; (72) 



in terms of the normalized rotation parameter, f N . 
3.5 Bifurcation points 

Given a homeidally striated ellipsoid, it has been shown above (see also C06) 
that the adjoint configuration, characterized by isotropic peculiar velocity 
distribution and normalized rotation parameter, (i>N)iso> exhibits same shape 
as a classical Jacobi ellipsoid of equal mass and axes, in real or imaginary 
rotation, where v = (fN)iso- Accordingly, the bifurcation point from ax- 
isymmetric to triaxial configurations is related to an axis ratio, £31, which is 
independent of the amount of systematic rotation and residual motion ex- 
cess, conform to Eq. (53). To gain more insight, let us equalize the alternative 
expressions of Eq. (46). The result is: 

M - ^44, = A 2 - C -^ej 2 A 3 ; (73) 

S>33 S33 

and the combination of Eqs. (59) and (73) yields: 



A 2 -A 1 e 



, , -f^3 ; (74) 

1 — e 21 e 21 

then the bifurcation points occur at a configuration, where the axis ratio, 
e 31 , is the solution of the transcendental equation (Caimmi 1996a): 

iim 4i^r T~ = e 3i^3 ; (75) 

£21— >i 1 — e 2 i 

where Eqs. (74) and (75) coincide with their counterparts related to isotropic 
residual velocity distributions (Caimmi 1996a,b). 

The contradiction with earlier results (Wiegandt 1982a,b; CM05) is ex- 
plained in the following way. Let us suppose that the generalized anisotropy 
parameters, Cn?C22,C33? can De arbitrarily fixed regardless from Eqs. (59), 
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and assume Cn = C22 also for triaxial configurations. Accordingly, Eq. (73) 
reads: 

A A f _ 

-31 



^£ = £4*; (76) 



1 — e 21 C33 e 21 

and the bifurcation points occur at a configuration where the axis ratio, 631, 
is the solution of the transcendental equation: 

2 A 2 — Ai (11 2 

hm e 21 - — = — e 31 A 3 ; (77) 

Qi-i 1 - e 21 ^33 

which coincides with Wiegandt (1982b) criterion for bifurcation, with re- 
gard to homeoidally striated Jacobi ellipsoids in rigid rotation. For a formal 
demonstration, see Appendix B. Then Wiegandt's criterion for bifurcation, 
expressed by Eq. (77), is in contradiction with Eqs. (59), contrary to the cur- 
rent one, expressed by Eq. (75). 

4 Transitions from and towards homeoidally 
striated Jacobi ellipsoids 

In a classical paper, Thuan & Gott (1975) modelled virialized elliptical galax- 
ies and related progenitor density perturbations at turn around as MacLaurin 
spheroids and rigidly rotating homogeneous spheres, respectively, and found 
a relation between the initial ratio of rotation to potential energy and the 
final shape. Surprisingly, the model succeeded in explaining why elliptical 
galaxies more flattened than E7 cannot exist. 

In a recent attempt (CM05) the above mentioned procedure has been 
generalized passing from MacLaurin spheroids to homeoidaly striated Jacobi 
ellipsoids, and the interested reader is addressed therein for deeper insight. 
What is relevant for the current investigation, shall be mentioned and further 
developed here. 

Let the initial and the final configuration be homeoidaly striated Jacobi 
ellipsoids, and let the former be denoted by a prime. The identity: 

rp _ -^rot -t^rot -^scl / 771 \ . 

TOt Z7V 771/ 771 V ^SeU ' 

A-ot --^sel -^sel 
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owing to Eqs. (5) and (18a), may be cast under the equivalent form (CM05): 
a[ _ (3 3 M S 11* £ Iot (7ga) 

(78b) 



ai ft S> U £U 
J ■ ft - M 



E, 



rot 



E s , 



involving dimensionless quantities only. 
On the other hand, the identity: 

E = E sel + E rot + Ej 
E 



E 
~E' 



E' 



^( E Ll + E TOt + E ' P CC + K SC ) 5 

by use of Eqs. (23), and (26), may be cast under the equivalent form (CM05): 



K 



sel 



E Ll 



E. 



sel 



rot 



E 'scl 



1 E se \ 1 i? se i c 



Pi 



— _ £ TO t — £pcc ~~ £qsc) j 

E „ E rot 



E> 
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sel 
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E 



sel 



f — 



c = 



E n 



E se \ 
Esel + 2E rot 



2(E osc + E ] 



(79a) 

(79b) 
(79c) 



pec j 



where E osc , E pec , represent the kinetic energy of systematic radial and non 
systematic motions, respectively. 

The combination of Eqs. (5), (18a), and (79) yields after some algebra: 



n2 S a 'l 



2C - 1 C - i 



2C 



c 



£ 



rot 



-/9e(1 £ ro t £pec ^osc) ~~ 5 



(80) 



involving dimensionless quantities only. 

The substitution of Eq. (78a) into (80) produces a second-degree equation 
in £ r ' ot , as (CM05): 



f' 2 - 2hf 

'-'rot z,L, °rot 



c = ; 



(81a) 
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M 



S\ 2 TZ' 

— I -^c-rot 



/?|/5e \S'J K 



2C-1 C - 1 



2C 



c 



rot 



(81b) 
(81c) 



the (reduced) discriminant of this equation is (CM05): 



b 2 - 



P 



M 



~ I "^-C-rot 



ffifc \S'J K 



2C - 1 C - i 



2C 



c 



rot 



(82) 



with regard to a transition from an initial to a final configuration, where all 
the parameters which appear in Eq. (82) are specified, except the axis ratios, 
e 2 i and 631. The condition, A = 0, via Eq. (82), represents a curve in the 
(Oe 2 ie3i) plane. The transition is forbidden for all values of the axis ratios, 
which make a negative discriminant, and then imaginary solutions. 
The solutions of Eq. (81a) are: 



3ut = &=F W 



Pm 



s\ 2 n' 



x 



2C - 1 C - 1 



£ 



rot 



1/2 



2C C 

and the combination of Eqs. (78a) and (83) yields: 



(83) 
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Pm & ^ro 



X-i C-i 



2C 
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b 2 - 



-V— £ 
ffifc \S'J K 

1/2 



rot 



(84) 



where the rotation parameter, £ r ' ot , and the axis ratio, ai/a[, are left un- 
changed for different departures from mass, angular momentum, and energy 
conservation, provided the ratios: 
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Pf 
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PIP* 
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p] 

Pm 



(85) 



do not vary. 
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5 The ending point of the sequence of ellip- 
tical galaxies 

The results of the current paper allow an extension and a generalization of the 
classical physical interpretation of the Hubble (1926) sequence (e.g., Jeans 
1929, Chap. XIII, §§ 298-303) to anisotropic peculiar velocity distributions. 
Owing to the results of Sect. 3, it will suffice to restrict to solid-body ro- 
tating configurations with isotropic peculiar velocity distributions, but with 
imaginary rotation also taken into consideration. 

5.1 Classification of galaxies 

It is worth recalling that, before establishing their similarity to the Milky 
Way, galaxies were referred to as "great nebulae" or "nebulae" . In the words 
of Jeans (1929, Chap. XIII, §298): 

"Hubble finds that it is not possible to place all observed nebulae 
in one continuous sequence; their proper representation demands a 
Y-shaped diagram. (...) 

The lower half of the Y is formed by nebulae of approximately 
elliptical or circular shape. These are subdivided into eight classes, 
designated E0, El, E7, the numerical index being the integer 
nearest to 10[(a — b)/a\, where a and b are the greatest and the least 
diameter of the nebulae as projected on the sky. Thus E0 consists 
of nearly circular nebuale (b > 0.95a), while E7 consists of nebulae 
for which b is about 0.3a, this being the greatest inequality of axes 
observed in the "elliptical" nebulae. (...) 

The upper half of the Y-shaped diagram consists of two distinct 
branches, one of which is found to contain a far larger number of 
nebulae than in the other. The principal branch contains the nor- 
mal "spiral" nebulae, which are characterized by a circular nucleus 
from which emerge two (or occasionally more) arms of approximately 
spiral shape. (...) These nebulae are subdivided into three classes, 
designated Sa, Sb, Sc, class Sa fitting almost continuously on to 
class El. 

The minor branch contains a special class of spirals characterized 
by the circumstance that the spiral arms appear to emerge from the 
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two ends of a straight bar-shaped or spindle-shaped mass. (...) 

About 97 per cent of known extra-galactic nebulae are found to 
fit into this Y-shaped classification. The remaining 3 per cent, are 
of irregular shape, and refuse to fit into the classification at all. (...) 
The irregular nebulae are distinguished by a complete absence of 
symmetry of figure and also by the absence of any central nucleus. 
(...) 

Apart from the irregular nebulae, Hubble states that, out of more 
than a thousand nebulae examined, less than a dozen refused to fit 
into the Y-shaped diagram at all, while in less than 10 per cent of 
the cases was there any considerable doubt as to a proper position 
of a nebula in the diagram. Clearly, the Y-shaped diagram provides 
a highly satisfactory working classification." 

A more detailed description of the classification scheme was performed 
later (Hubble 1936), and subsequently modified as the collection of large- 
scale plates of giant galaxies grew. Reproductions of many of these plates 
were published posthumously (Sandage 1961) and the introduction therein 
is generally regarded as the definitive exposition of the Hubble scheme. At 
the junction of elliptical and spiral galaxies, comes a class of galaxies (not 
mentioned in early classifications) known as lenticulars. These galaxies are 
designated as type SO or type SBO according to whether or not they are 
barred. For further details refer to e.g., Mihalas & Binney (1981, Chap. 5, 
§5.1). 

The Hubble system, as conceived by Hubble and further developed by 
Sandage, is defined in terms of type-examples which are almost exclusively 
giant galaxies (M pg < —19). However, the most numerous type of galaxy in 
the Local Group is either dwarf elliptical or dwarf spheroidal. The surface 
brightness of a typical object belonging to the former class, is perfectly nor- 
mal for an elliptical galaxy, so that it differs from a giant elliptical galaxy 
only in linear size and absolute magnitude. Dwarf spheroidal galaxies, on the 
other hand, are very low surface brightness objects. Thus it may safely be 
expected that dwarf galaxies are the most numerous objects in the universe. 

Hubble's original scheme has generally been considered satisfactory in 
regard to the ellipticals, but it has been said that Hubble's classification of 
the spirals is incomplete and that his treatment of irregular galaxies was 
quite inadequate. A few of the more important attempts at reclassification 
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of types later than E (that is, to the right side of the ellipticals in Hubble's 
diagram) are described below. 

Hubble's two-dimensional scheme has been made three-dimensional to 
include explicit reference to rings and s-shaped objects, and, in addition, the 
sequence has been extended to Sd-m and Im galaxies (de Vaucouleurs 1959). 
The Sd class overlaps Hubble's Sc class to some extent, but it also contains 
more extreme objects which are classified as Type I Irregulars (Irr I) in 
Hubble's scheme. The Sm and Im classes contain the remaining galaxies of 
Hubble's Irr I class. 

An alternative scheme (the Yerkes system) relies on a fundamental pa- 
rameter designating the population group or concentration class of a galaxy 
(McClure & van den Bergh 1968; Morgan et al. 1975). This parameter runs 
from k to a such that, among normal (in the sense of "non active") galaxies, 
those of type k have the highest degree of central concentration of their light, 
and those of type a have the smallest central bulges and the most diffuse light 
distributions. Galaxies of type gk, g, fg, f, af, have various intermediate 
light concentration. 

Many features of both the classical Hubble system and the Yerkes system 
are incorporated in a new scheme (the DDO system), where the lenticular 
galaxies are arranged parallel to the spirals rather than before them, and a 
new class of galaxy, the "anemic" spiral, is interposed between the spirals 
and the lenticulars (van den Bergh 1960a,b; 1976). The resulting figure is a 
trident diagram, and the related sequences of barred galaxies are described 
by a separate trident. 

For further details and references on the classical Hubble system, the 
de Vaucouleurs system, the Yerkes system, and the DDO system, see e.g., 
Mihalas & Binney (1981, Chap. 5, §5.1). A deep analysis on the physical 
parameters along the Hubble sequence may be found in e.g., Roberts & 
Haynes (1994). The variation of gas content along the Hubble sequence may 
be explored using a new catalogue of normal (in the sense of "isolated") 
galaxies (Bettoni et al. 2003). 

5.2 Physical interpretation 

Soon after the early Hubble (1926) classification of the "great nebulae", a 
physical interpretation was provided by Jeans (1929, Chap. XIII, §299): 
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Obviously the proper physical interpretation of the classification 
just described is of the utmost importance to cosmogony. 

A first and most important clue is provided by the fact that 
numbers of the great nebulae are known to be in rotation. (...) 

The symmetry of figure shewn by nebulae of the E and Sa types 
is precisely such as rotation might be expected to produce, and this 
suggests an inquiry as to how far the observed figures of the regular 
nebulae can be explained as the figures assumed by masses rotating 
under their own gravitation" . 

The conclusion is (Jeans 1929, Chap. XIII, §302): 

"Remembering that rotation has actually been observed in a num- 
ber of nebulae, there seem to be strong reasons for conjecturing that 
the observed configuration of the nebulae may be explained in general 
terms as the configurations of rotating masses" . 

More specifically, two classes of models are considered, where density 
profiles range between the limiting cases (i) homogeneous configurations i.e. 
classical MacLaurin spheroids and Jacobi ellipsoids, and (ii) mass points sur- 
rounded by massless atmospheres i.e. Roche systems. Rigidly rotating mass 
distributions of the kind considered may attain any shape between the ex- 
treme boundaries (a) spherical i.e. nonrotating, and (b) centrifugal support 
at the end of major equatorial axis. With respect to the latter configuration, 
further rotation implies equatorial shedding or top major equatorial axis 
streaming of matter, for axisymmetric and triaxial configurations, respec- 
tively. For a fixed density profile, the sequence of rigidly rotating configura- 
tions starts from the spherical shape and ends at the shape where centrifugal 
support is first attained. 

The above interpretation suffers from two main points (Jeans 1929, Chap. XIII, 
§§300-302). First, the most flattened elliptical configuration (if related to a 
figure of revolution), should occur when centrifugal support takes place at 
the bifurcation point, from axisymmetric to triaxial shapes. On the other 
hand, the above mentioned limiting configuration is less flattened than E7. 
Second, triaxial bodies of the kind considered cannot be figures of equilib- 
rium if a great central condensation of mass is present. On the other hand, 
it is the case for barred spirals. 
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Owing to the results of Sect. 3, homeoidally striated Jacobi ellipsoids 
with arbitrary peculiar velocity field, and systematic motions reduced to ro- 
tation around a fixed axis, may be related to classical Jacobi ellipsoids, pro- 
vided imaginary rotation is taken into consideration. Then the occurrence 
of anisotropic peculiar velocity distributions does not affect the physical in- 
terpretation under discussion, and the questions to be clarified remain the 
above mentioned two. 

Accordingly, the second sentence quoted from Jeans at the beginning of 
the current Section, could be generalized as: "A first and most important 
clue is provided by the fact that galaxies are known to be in (real or imagi- 
nary) rotation". In particular (i) spiral galaxies apparently rotate; (ii) faint 
(Mb > —20.5) elliptical galaxies (Davies et al. 1983), SA bulges (Kor- 
mendy & Illingworth 1982), SB bulges (Kormendy 1982), and a few bright 
(Mb < —20.5) elliptical galaxies (Illingworth 1981) appear to be supported 
by systematic rotation; (iii) the majority of bright elliptical galaxies (Illing- 
worth 1981) appear to be supported by anisotropic peculiar velocity distri- 
bution; and (iv) elliptical-like galaxies with dust lanes (Sharpies et al. 1983), 
similarly to bright ellipticals, appear to be supported by either systematic 
rotation or anisotropic peculiar velocity distribution. For a more detailed 
discussion refer to Caimmi (1983). 

In absence of a unified theory of systematic and random motions, with 
regard to rotation around a fixed axis, centrifugal support and anisotropic 
pressure must necessarily be thought of as independent contributions to the 
shape of the system (e.g., Binney 1976, 1978, 1980). On the contrary, the 
current attempt comes back to Jean's conception, that the shape of galaxies 
is determined by rotation, conceived as real (related to centrigugal support) 
or imaginary (related to anisotropic peculiar velocity distribution). To this 
aim, our attention shall be restricted to the Hubble classification for the 
following reasons. First, it is defined in terms of almost exclusively bright 
galaxies, where both centrifugal support and anisotropic peculiar velocity 
distribution have been observed. Second, it has generally been considered 
satisfactory in regard to ellipticals. 

The detection of highly flattened ellipticals is expected to be rare, as they 
must necessarily be viewed edge-on. In addition, the advent of more refined 
instruments and techniques has led to different classifications. For instance, 
NGC 3115 is quoted as E7 in Jeans (1929, Chap. XIII, §298, Plate IX), 
E7 / SO in Sandage's Hubble Atlas of Galaxies (e.g., Mihalas & Binney 1981, 
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Chap. 5, § 5.1, Fig. 5.3), and SO in Larsen et al. (1983); NGC 3377 is quoted 
as EG in Sandage's Hubble Atlas of Galaxies (e.g., Mihalas & Binney 1981, 
Chap. 5, §5.1, Fig. 5.3), E5-6 in Copin et al. (2004); and the junction of the 
elliptical and spiral galaxies occurs at E7 type in the Hubble's fork diagram 
(e.g., Mihalas & Binney 1981, Chap. 5, §5.1, Fig. 5.2) and at E6 type in 
the van den Bergh's trident diagram (e.g., Mihalas & Binney 1981, Chap. 5, 
§5.1, Fig. 5.8). Additional examples of highly flattened (E6) ellipticals are 
NGC 821 and NGC 4697 which, together with NGC 3377, are known to host 
supermassive (M « 1O 8 M ) black holes (Soriaet al. 2006). Different features 
exhibited by highly flattened elliptical and lenticular galaxies, are shown by 
typical objects belonging to each class, as in the cases shortly reported below. 

NGC 3377 is a prototypical "disky" elliptical galaxy with "boxy" outer 
isophotes. It has a power-law central luminosity profile and its total abso- 
lute magnitude of about —19(5) is intermediate between that of the classical 
"boxy" giant ellipticals and "disky" lower-luminosity objects. Both dynam- 
ical model and the M* — a relation suggest the presence of a massive black 
hole. For further details and references refer to e.g., Copin et al. (2004); 
Soria et al. (2006). 

NGC 3115 has long been assumed to be the prototype of the SO galaxy 
type: a bulge-dominated galaxy with an embedded disk and very little gas 
and dust. The system has a nearly edge-on inclination and contains a double 
disk structure with an outer Freeman type II disk, which exhibits a weak 
spiral arm structure, and a nuclear disk of size about fifty times shorter. The 
spheroidal component does not seem to follow the classical r 1 / 4 law, and the 
flattened halo extends up to nine times outside the larger disk. The bulge 
appears to be supported by systematic rotation and evidence has been found 
for the presence of both a central dark mass and a massive dark halo. For 
further details and references refer to e.g., Emsellem et al. (1999). 

5.3 The E sequence within Ellipsoidland 

In dealing with a physical interpretation of the early Hubble sequence, it is 
convenient to define axis ratios of intrinsic configurations in a different way. 
Given a homeoidally striated Jacobi ellipsoid, let a, b, c, be the semiaxes, 
where a > b > c without loss of generality. Accordingly, e ca < < 1; in 
addition, the minor and the major axis coincide with the rotation axis for 
oblate-like and prolate-like configurations, respectively. 
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The whole range of possible configurations in the (Oeft a e ca ) plane defines 
Ellipsoidland (the term is from Hunter & de Zeeuw 1997). Ellipsoidland 
is a triangle where two sides are of unit length and an angle is right, as 
shown in Fig. 1. Oblate configurations lie on the vertical side, €b a = 1, the 
spherical configuration (square) is on the top, e ca = 1, and the flat circular 
configuration is on the bottom, e ca = 0. Prolate configurations lie on the 
inclined side, e ca = e& a , the spherical configuration is on the top, e ca = 
€b a — 1, and the prolate oblong configuration is on the bottom, 

tea Cfca 

0. Flat configurations lie on the horizontal side, e ca = 0, the flat circular 
configuration is on the right, tb a = 1, and the flat oblong configuration is 
on the left, tb a = 0. It is worth noticing that flat-oblong and prolate-oblong 
configurations are not coincident (e.g., Caimmi 1993; CM05). Non flat, non 
axisymmetric configurations, lie within Ellipsoidland. 

The axis ratio correlation, e ca vs. e^, related to adjoint configurations to 
homeoidally striated Jacobi ellipsoids, is represented within Ellipsoidland in 
Fig. 1. The bifurcation point from MacLaurin spheroids to Jacobi ellipsoids 
is marked by a triangle. The bifurcation point from MacLaurin spheroids and 
Jacobi ellipsoids to pear-shaped configurations is marked by a St. Andrew's 
cross and a Greek cross, respectively. With regard to real rotation, the cor- 
relation is the known one involving classical MacLaurin spheroids and Jacobi 
ellipsoids. With regard to imaginary rotation, the correlation reads e ca = o, a , 
owing to lack of bifurcation points. Then configurations in real and imagi- 
nary rotation branch off from the nonrotating spherical configuration. The 
former sequence proceeds down along the oblate side of Ellipsoidland, until 
the bifurcation point is attained and the curve enters Ellipsoidland, finishing 
when the next bifurcation point is also attained. The latter sequence pro- 
ceeds down along the prolate side of Ellipsoidland, until the prolate-oblong 
configuration is attained, never entering Ellipsoidland. On the other hand, 
a sequence may stop earlier, when centrifugal support at the ends of major 
equatorial axis is attained, which depends on the density profile. In fact, it 
makes the sole difference between homeoidally striated Jacobi ellipsoids and 
their adjoint counterparts (for further details, see CM05). 

With regard to Ellipsoidland, the observed lack of elliptical galaxies more 
flattened than E7 translates into the inequality, e ca > 0.25, the threshold 
being represented by the lower dashed horizontal line in Fig. 1. The locus of 
edge-on configurations, viewed along a direction coinciding with equatorial 
minor axis, where 0.25 < e ca < 0.35, and then belonging to class E7, is 



30 



1.5 



1.0 
» 0.5 

0.0 

-0.5 

-0.5 0.0 0.5 1.0 1.5 

Figure 1: Axis ratio correlation within Ellipsoidland, with regard to adjoint 
configurations to homeoidally striated Jacobi ellipsoids. The semiaxes are a, 
b, c, where a > b > c or e ca <J ei, a <J 1, without loss of generality. The locus of 
flat configurations is e ca — 0, < e^a < 1- The locus of oblate configurations 
is €ba = 1, < e ca < 1. The spherical configuration (square) is defined as 
e ca = e& a = 1. The oblong configuration is defined as e ca = €b a = 0. The locus 
of prolate configurations is e ca = e& a , < e^a < 1- Different configurations 
occur in the flat oblong limit, e ca = 0, e^a — ► 0, and in the prolate oblong 
limit, e ca — >• 0, efe a — > 0, e ca = efe a . A band, bounded by two horizontal dashed 
lines, defines class E7 in Hubble (1926) classification i.e. 0.25 < e ca < 0.35 
for edge-on ellipsoids where the line of sight coincides with the direction of 
the minor equatorial axis. The bifurcation point from MacLaurin spheroids 
to Jacobi ellipsoids is marked by a iangle. The bifurcation point from 
MacLaurin spheroids and Jacobi ellipsoids to pear-shaped configurations is 
marked by a St. Andrew's cross and a Greek cross, respectively. 




represented in Fig. 1 as a band bounded by two dashed lines, e ca = 0.25 
and e ca = 0.35, respectively. A change in direction of the line of sight could 
project an intrinsic ellipsoid within the above mentioned band, on a class Ei, 
i < 7. 

An inspection of Fig. 1 shows that the bifurcation points from triaxial and 
axisymmetric configurations, towards pear-shaped configurations, are consis- 
tent, or marginally consistent, with the E7 band on Ellipsoidland. Then the 
occurrence of triaxiality provides a physical interpretation to the observed 
absence of elliptical galaxies (and spiral bulges) more flattened than E7, 
with regard to the oblate-like branch of the sequence. On the other hand, 
the above interpretation cannot apply to the prolate branch of the sequence, 
as no bifurcation point has been found therein, and some kind of instability 
must be considered. 

The absence of elliptical galaxies more flattened or elongated than E7 
might be due to bending instabilities, as suggested long time ago from an- 
alytical considerations involving homogeneous (oblate and prolate) spheroids 
(Polyachenko & Shukhman 1979; Fridman & Polyachenko 1984, Vol. 1, Chap. 4, 
Sect. 3.3, see also pp. 313-322; Vol. 2, p. 159) and numerical simulations in- 
volving inhomogeneous (oblate and prolate) spheroids (Merritt & Hernquist 
1991; Merritt & Sellwood 1994). The amount of figure rotation seems to be 
unimportant to this respect (Raha et al. 1991; Merritt & Sellwood 1994). 

This conclusion is supported by recent results from iV-body numerical 
simulations, where spherically symmetric, unstable, radially anisotropic, one- 
component 7 models were taken as starting configurations (Nipoti et al. 
2002). The related end-products, in accordance with previous results (e.g., 
Merritt & Aguilar 1985; Stiavelli & Sparke 1991), were found to be in general 
prolate systems less flattened than E7 (Nipoti et al. 2002, Fig. 1 therein). 

The theoretical explanation of the result, in terms of a dynamical bending 
instability (Merritt & Sellwood 1994), is generally recognized to explain also 
the maximum elongation of simulated nonbaryonic dark matter haloes (e.g., 
Bett et al. 2007). 

5.4 Tidal effects from hosting dark matter haloes 

Current ACDM cosmologies, which provide a satisfactory fit to data from 
primordial nucleosynthesis and cosmic background radiation, predict galaxies 
are embedded within dark matter haloes. Then it is a natural question 
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to what extent the presence of hosting dark matter haloes may affect the 
above interpretation of the early Hubble sequence. To this aim, an idealized 
situation shall be analysed. 

Let us represent elliptical galaxies and their hosting dark haloes as con- 
centric and coaxial classical Jacobi ellipsoids, one completely lying within 
the other. Accordingly, the two bodies must necessarily rotate at the same 
extent, and/or one of them (the outer in the case under discussion) has to 
be axisymmetric. It can be seen that the effect of the outer ellipsoid on the 
inner one is to shift bifurcation points towards more flattened configurations 
with respect to a massless embedding subsystem (Durisen 1978; Pacheco et 
al. 1986; Caimmi 1996a). 

More specifically, the axis ratio of the configuration at the bifurcation 
point is the solution of the transcendental equation (Caimmi 1996a): 



[(^)3l] 2 (A) 



i)3 



4[(e,) 31 ] 2 , 4-4[(e 4 ) 3 i 12 



+ K 



3-2[(e,) 31 ] 2 3-2[( Ci ) 3 i 



(86a) 



Mi (aj fc ymy3 pi 

k = y-^ ; < k < +oo ; 86c 

1/11/21/3 (A)3 

where the indices, i, j, denote inner and outer ellipsoid, respectively. In the 
special case of massless outer ellipsoid, k — 0, Eq. (86a) reduces to (75). 

Further analysis shows that the embedded configuration, at the bifurca- 
tion point, can be as flattened as E7 provided the parameter, k, is close to 
unity. On the other hand, the related parameters are to be consistent with 
observations and cosmological models. To this aim, inhomogeneous density 
profiles must be considered. 

The generalization of Eq. (86a) to homeoidally striated Jacobi ellipsoids 
demands to restart from a generalized formulation of the virial theorem, 
where the tidal potential is also included (Brosche et al. 1983; Caimmi et al. 
1984; Caimmi & Secco 1992). The repetition of the same procedure used in 
the current paper, yields: 

- = — (87) 

2/12/22/3 (^)sel {Ai) 3 



33 



where (^)tid is an additional factor in the expression of the potential tidal 
energy, related to the tidal action of the embedding subsystem on the em- 
bedded one (Caimmi 2003; CM05). 

It may safely be expected that the external boundary is less flattened than 
the internal one and more flattened than a sphere. Accordingly, the shape 
factor ratio appearing in Eqs. (87), ranges as 1/2 < (Aj) 3 /(Ai) 3 < 1 for oblate 
configurations, and 1 < (A,-) 3 / (Aj) 3 < 10/3 for prolate configurations. Then, 
to a first extent, (A,) 3 / (Ai) 3 pa 1, i.e. the boundaries are similar and similarly 
placed ellipsoids. With this restriction, the factor, (i^-)tid, is a profile factor, 
which may be expressed by a simple formula (Caimmi 2003): 

Mt* = -l - ( , m f\ ™ {ext) - ; (88) 

O I v i ) mas I Vj ) mas Ho 

where y is the scaling radius ratio of outer to inner subsystem, and w^'^ is 
an additional profile factor. 

Typical elliptical galaxies and their hosting dark haloes may safely be 
represented as homeoidally striated Jacobi ellipsoids, where the star and 
dark subsystem are described by generalized power-law density profiles of 
the kind: 

fu(Q = <a v I X = — — I u = i,j ; (89a) 

f = r " • » = Ru ■ (89b) 
{r )u {ro)u 

& = Votj ; y^i = y ^j ; (89c) 

s = f : (89d) 

according to Eqs. (1), and the choices (a, (3, 7) = (1,4,1) (Hernquist 1990) 
and (a, /3, 7) = (1,3,1) (Navarro et al. 1995, 1996, 1997) are adopted to 
describe the star and the dark subsystem, respectively. The values of input 
and output parameters of the model are listed in Tab. 1, with regard to a 
ACDM cosmology where Vt M = 0.3, fi A = 0.7, tt b = 0.0125/?r 2 , h = 2^ 2 . 
For further details and references, see CM03, CM05. 

The presence of a massive halo has little influence on the location of the 
bifurcation point from axisymmetric to triaxial configurations, which is found 
to occur at an axis ratio, (e 3 i)bif = 0.580088, related to k — 0.00441624. 
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value 


output 


value 




40/3 


( Vi ) mas 


10.38399 




10 


( ) mas 


17.86565 


(r )j/kpc 


3.21 


(^)sel 


1.44444 


(r ) i /kpc 


36.20 


(^i)scl 


0.6268271 


Ri/kpc 


42.77 


(^j)tid 


0.3518317 


Rj /kpc 


362.04 


(^i)tid 


0.2217760 


Mi/1O 1O M 


8.33 


(^ijrot 


1/3 


Mj/10 10 M Q 


91.67 


(^jjrot 


1/3 




11.29 


{Vi)rot 


1/2 


y 


8.47 


(Vjjiot 


1/2 


m 


11 


w (est) 
w (int) 

( ^ij )int 

[yjijint 
( 2r) rot v rot \ 
\ u scl J i 
( 2r?rotfrot A 
V ^scl / j 
m (Vij)tid 
V 3 (f»)sel 
(C3l)bif 


-0.3955800 
-3.564028 
1.760852 
1.760852 
0.230769 
0.531779 

0.00441624 
0.580088 



Table 1: Values of input and output parameters in modelling elliptical galax- 
ies and their hosting dark haloes as similar and similarly placed, homeoidally 
striated Jacobi ellipsoids, characterized by Hernquist (1990) and Navarro 
et al. (1995, 1996, 1997) density profiles, respectively. The parameters of 
the related ACDM cosmology have been chosen as Vt M = 0.3, Vl\ = 0.7, 
VL b = 0.0125/i~ 2 , h = 2 -1 / 2 . For further details on the output parameters see 
e.g., CM03. 
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A similar result is expected to occur for the bifurcation point from both 
axisymmetric and triaxial to pear-shaped configurations. Then the above 
mentioned points continue to be consistent, or marginally consistent, with 
the El band on Ellipsoidland, even in presence of (typical) massive dark 
haloes, concerning the oblate-like branch of the sequence. On the contrary, 
it is suggested the existence of some kind of instability, which does not allow 
prolate configurations in rigid imaginary rotation, more elongated than E7, 
even in presence of (typical) massive dark haloes, according to recent results 
from iV-body simulations (Nipoti el al. 2002). 

5.5 Cosmological effects after decoupling 

About thirthy years ago, Thuan & Gott (1975) idealized elliptical galax- 
ies as MacLaurin spheroids, resulting from virialization after cosmological 
expansion and subsequent collapse and relaxation of their parent density 
perturbation. A generalization of the method to triaxial configurations and 
anisotropic peculiar velocity distributions, has been performed in CM05, and 
an interested reader is addressed therein for further details. Owing to the 
results of Sect. 3, the case of isotropic peculiar velocity distribution, involving 
both real and imaginary rotation, can be considered without loss of general- 
ity. 

Our attention shall be limited to dark matter haloes hosting giant galax- 
ies, as massive as about 1O 12 M . Accordingly, it is assumed 5 rec = 0.015 
as a typical overdensity index of the initial configuration, taken to be at 
recombination epoch, and a relaxed final configuration i.e. ( = 1. 

With regard to the initial configuration, the following approximations 
hold to a good extent: (i) spherical shape; (ii) homogeneous mass distri- 
bution; (iii) negligible rotation energy; (iv) negligible peculiar energy. The 
changes in mass, angular momentum, and total energy, respectively, during 
the transition from the initial to final configuration, are expressed by the 
parameters, (3m, Aj, /3e, defined by Eqs. (78b) and (79b), where the initial 
configuration is marked by the prime. For further details, see CM05. 

Though negligible with respect to the potential and expansion energy, 
the rotation energy of the initial configuration affects the shape of the fi- 
nal configuration, as described by Eqs. (78), (81), (83), (85), where different 
changes in mass, angular momentum, and tidal energy, are related to a same 
final shape for an assigned initial configuration, provided the parameter, fls , 
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defined by Eq. (85), does not vary. The choice, Ps = 1/3600, in particular 
{Pm, Pj, Pe) = (1,60, 1), holds for dark matter haloes hosting giant elliptical 
galaxies, and it shall be assumed here. For further details, see CM05. 

The axis ratios of the final configuration, €31 and function of the 

parameter, ks = £' TO J{\ — £' osc — £L C ), are plotted in Fig. 2, with regard to 
a Navarro et al. (1995,1996,1997) density profile (S = 9.20678; see CM05, 
Tab. 1 therein, for further details) in rigid rotation. Each curve is symmetric 
with regard to a vertical axis, ks = 0.5, where a local minimum is attained. 
The plane, (Orege), may be divided into three regions, namely (a) 0.163190 < 
Ks < 0.836810, where triaxial configurations occur, ranging from (€21,631) = 
(1,0.582724) to (e 2 i,e 3 i) = (0.363174,0.302297), the latter related to k £ = 
0.5; (b) < k 8 < 0.163190, 0.836810 < k £ < 1, where oblate configurations 
occur, ranging from (e 2 i,e 3 i) = (1,0.582724) to (e 2 i, e 31 ) = (1,1), the latter 
related to Ks — (bound), k £ = 1 (unbound), respectively; (c) —00 < 
K s < 0, 1 < Ks < +00, where prolate configurations occur, ranging from 
(e 2 i,63i) = (1,1) to (e 2 i,e3i) = (l,+oo), the latter related to k £ — > —00 
(bound), k £ — > +00 (unbound), respectively. It is worth recalling that oblate- 
like (e 3 i < e 2i < 1) and prolate-like (e 2 i < 1 < e 3 i) configurations, are related 
to real and imaginary rotation, respectively. 

An inspection of Fig. 2 shows that, in the case under discussion, cosmolog- 
ical effects due to expansion prevent dark matter haloes in real rotation from 
being more flattened than E7 (dashed horizontal band). The same holds for 
embedded giant elliptical galaxies, provided the related shapes may safely be 
thought of as similar and similarly placed. On the other hand, an arbitrary 
flattening can be attained by dark matter haloes in imaginary rotation, un- 
less some kind of instability occurs, which similarly prevents configurations 
more elongated than El. 

As outlined in Subsect. 5.3, bending instabilities have been suggested long 
time ago as a viable mechanism to this respect (Polyachenko & Shukhman 
1979; Fridman & Polyachenko 1984, Vol. 1, Chap. 4, Sect. 3.3, see also pp. 313- 
322; Vol.2, p. 159; Merritt & Hernquist 1991; Merritt & Sellwood 1994). 

The above results depend on both density profile and rotation velocity 
profile of the final configuration, via the coefficient, c, appearing in Eq. (81c). 
Owing to Eqs. (7) and (18c), Eq. (81c) translates into: 

ram 
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Figure 2: Equatorial (upper curve) and meridional (lower curve) axis ratio 
of the final configuration, as a function of the parameter, ks = £' TO J{\ — 
£' osc — £p CC ), for sequences of relaxed Navarro et al. (1995, 1996, 1997) 
density profiles in rigid rotation, related to dark matter haloes hosting gi- 
ant elliptical galaxies. The curves are symmetric with respect to a vertical 
axis, K£ = 0.5, and a minimum extremum point occurs at (ks, €21, €31) = 
(0.5,0.363174,0.302297), corresponding to the most flattened, oblate-like 
configuration which is allowed. Values of k £ within or outside the range, 
< Ks < 1, are related to real or imaginary rotation, respectively. Values of 
Ks below or not below unity, are related to bound (finite extent) and unbound 
(infinite extent) configurations, respectively. Oblate-like (631 < e 2 ± < 1) 
triaxial configurations lie between the inner, vertical dotted lines. Oblate 
(f3i < £21 = 1) axisymmetric configurations lie between the inner and the 
corresponding outer, vertical dotted lines. Prolate (e 2 i = 1 < e 3 i) axisym- 
metric configurations lie outside the outer, vertical dotted lines. The equa- 
torial axis ratio is e = £21. The polar axis ratio is e = 631 for oblate and 
e = 613 for prolate configurations. The horizontal dashed lines define class 
E7 in Hubble (1926) classification i.e. 0.25 < e < 0.35 for edge-on ellipsoids 
where the line of sight coincides with the direction of the minor equatorial 




u 2 B 5 u 2 

fa=Pe^ = ik-? S!L ; (90b) 

^ram PjPE ^ram 

where different changes in mass, angular momentum, total energy, density 
profile, and rotation velocity profile, are related to the same final shape for 
an assigned final configuration, provided /%y does not vary. The choice, j3e> = 
6.007519-10~ 6 , in particular (f3 M , Pj, Pe, v se \, ^m) = (1, 60, 1, 0.610045, 17.20783), 
holds for Fig. 2. 



6 Conclusion 

Elliptical galaxies have been modelled as homeoidally striated Jacobi ellip- 
soids where the peculiar velocity distribution is anisotropic, or equivalently 
as their adjoint configurations i.e. classical Jacobi ellipsoids of equal mass 
and axes, in real or imaginary rotation. Reasons for the coincidence of bi- 
furcation points from axisymmetric to triaxial configurations in both the 
sequences (CM06), contrary to earlier findings (Wiegandt, 1982a,b; CM05), 
have been presented and discussed. The occurrence of centrifugal support at 
the ends of major equatorial axis, has been outlined. 

The existence of a lower limit to the flattening of elliptical galaxies has 
been investigated in dealing with a number of limiting situations. More 
specifically, (i) elliptical galaxies have been considered as isolated systems, 
and an allowed region within Ellipsoidland (Hunter & de Zeeuw 1997), related 
to the occurrence of bifurcation points from ellipsoidal to pear-shaped con- 
figurations, has been shown to be consistent with observations; (ii) elliptical 
galaxies have been considered as embedded within dark matter haloes and, 
under reasonable assumptions, it has been shown that tidal effects from host- 
ing haloes have little influence on the above mentioned results; (iii) dark mat- 
ter haloes and embedded elliptical galaxies, idealized as a single homeoidally 
striated Jacobi ellipsoid, have been considered in connection with the cos- 
mological transition from expansion to relaxation, by generalizing an ear- 
lier model (Thuan & Gott 1975), and the existence of a lower limit to the 
flattening of relaxed (oblate-like) configurations, has been established. On 
the other hand, no lower limit has been found to the elongation of relaxed 
(prolate-like) configurations, and the observed lack of elliptical galaxies more 
elongated than E7 has needed a different physical interpretation such as the 
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occurrence of bending instabilities (Polyachenko & Shukhman 1979; Merritt 

6 Hernquist 1991). 
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Appendix 

A Some properties of ellipsoid shape factors 

Ellipsoid shape factors, A ± , A 2 , A 3 , obey the following inequalities (Caimmi 
1996a): 

a" Ax > a 2 A 2 > a 3 A 3 ; a x > a 2 > a 3 ; n > 2 ; (91a) 
a\A x < a n 2 A 2 < a^A 3 ; a ± > a 2 > a 3 ; n < 1 ; (91b) 

where inequalities (91a) and (91b), the latter restricted to n < 0, come 
from analytical considerations involving the explicit expression of the shape 
factors (e.g., MacMillan 1930, Chap. II, §33; A x = a~ 2 , A 2 = (3~ 2 , A 3 = 7" 2 , 
therein), and inequality (91b), restricted to < n < 1, results from an 
obvious generalization of a proof by Pacheco et al. (1989). 
Using Eqs. (6) and (7) yields: 

¥=4t ; ?=1 ' 2; (92) 

°33 e 3g /1 3 

and, owing to inequality (91a): 

A q >t 2 3q A 3 ■ a q >a 3 ; 2 = 1,2 ; (93) 

which implies S qq > S 33 for oblate-like configurations (a q > a 3 ), and S qq < 
S 33 for prolate-like configurations (a q < a 3 ). 
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In the limit of axisymmetric configurations, e 2 i — 1, 631 = e, A\ = A 2 = a, 
A3 = 7, and the following relations hold (e.g., Caimmi 1991, 1993): 



lima = ; lim 7 = 2 

£^0 e— 



a n 
hm — = — 

e-»o e 2 



lim a 



lim 7 = lim (57) = 



lim 



7 — a 



6-1 1 - e 2 
da 1 

"d7 = 1 -e 2 
d7 1 

~de~ ~ 1 -e 2 



o. a 

(l + 2e 2 )--2e 
e 

;i + 26 2 )^-- 



v = (^N)iso = a - e 2 7 ; 



(94) 

lim (e 2 7) = +00 ; (95) 

(96) 

(97) 

(98) 
(99) 



and keeping in mind the general property (e.g., Chandrasekhar 1969, Chap. 3, 
§17): 

A 1 + A 2 + A 3 = 2 ; (100) 

it can be seen that the first derivative of the rotation parameter, dv / de, is 
null provided the transcendental equation: 

6e 2 

(101) 



a = 



l + 8e 2 ' 

is satisfied. One solution is found to exist, which is related to oblate configu- 
rations. It corresponds to an extremum point where the function attains its 
maximum value (e.g., Chandrasekhar 1969, Chap. 5, §32). 



B Wiegandt criterion for bifurcation 

Given a collisionless, self-gravitating fluid in rigid rotation, where no internal 
energy transport occurs and the residual velocity is constant on the boundary, 
an upper limit for the point of bifurcation is (Wiegandt 1982a,b): 

^ 2 hi = V 12 , 12 ■ (102) 

where V pg . rs is the "super-matrix" : 

dV r 



J^p(x 1 ,x 2 ,x 3 )x p -^ 1 d 3 S ; (103) 
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/ n n f P{ x l, x 2i x 3)( X P X p){ x s x ' s ) ,3 c / ln A 

Vrsix^xi, x s ) = Gj s [{xi _ x , iY + {x2 _ x , 2 y + {xz _ x ^ fl , d 5,(104) 

which is expressed in terms of a generalized potential, V rs , and the integra- 
tions are carried over the whole volume, S, of the system. 

In the special case of homeoidally striated Jacobi ellipsoids, Eq. (102) 
identifies the exact point of bifurcation (Wiegandt 1982a,b) and the following 
relations hold (Wiegandt 1982b): 

Vpqwq = a (Esei) PP ', (105) 

Tip 



a 2 q A pq = ^—^ ; (106) 

€ pq 

where the products, a 2 q A pqi are shape factors which depend on the axis ratios 
only, similarly to A p , and symmetry with respect to the indices holds, A pq = 
A qp . In addition, Eq. (106) has been deduced from Chandrasekhar (1969, 
Chap. 3, §21), Eq. (107) therein. 

Though Wiegandt (1982b) analysis is restricted to binomial density pro- 
files (e.g., Perek 1962; Chandrasekhar 1969, Chap. 3, §20; Caimmi 1993), 
still it may be generalized to any kind of cored density profiles, defined as: 

Pw(0 = Pcf(Zw) ; /(0) = 1 ; p c = p(0) ; (107a) 
6y = ^ ; 0<Zw<l ; E w = l ; (107b) 

where the scaling radius and the scaling density are chosen to be equal to the 
radius, R, and the central density, p c , respectively. The following relations: 

P = —pw ; f = ; (108) 

Pc 

allow conversion from Eqs. (1) to Eqs. (107) and vice versa. 

Using Eqs. (107), the mass, the inertia tensor, and the potential-energy 
tensor, take the equivalent expression: 

M = (v mas ) w M c ; (109a) 
M c = — p c a?e 2 ie 3 i ; (109b) 
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Ipq = 8p q e pl e gl alM(is inr ) w ; (110) 
GM 2 

(E se l) pq = — (^scl)vK (-Bsel)pg ; (HI) 

and the comparison with Eqs. (2), (3), (4), (8), yields: 

(^mas)vK = ~ ^^mas ] (^inr)vK = ^inr ] (v S el)w = ^scl ] (H2) 

Pc " 

finally, using Eqs. (6) and (109)-(112), the potential-energy tensor takes the 
equivalent form: 



(E'seOpg = -knGp c Ip q A p ; (113a) 

4 Po ^scl^n 
3p c S 3 //,, 



^^^Wmas (n3b) 

> ~ 3 7 ■ 



being /c, by definition, a profile factor. 

Owing to Eqs. (8), (19), (21), (31), (32), (113b), the normalized rotation 
parameter, defined by Eq. (45), may be expressed as: 

l * = ^ ; (114) 

which coincides with the parameter, 0, defined in Wiegandt (1982b) in the 
special case of rigid rotation (z/ mr = z/ rot ; i] rot = 1/2; see CM05), as shown in 
Caimmi (1996b). Accordingly, Eq. (46) coincides with its Wiegandt (1982b) 
counterpart, Eq. (52) therein. 

The combination of Eqs. (102), (105), (113), (114), yields: 

v N = A 1 - a 2 2 A 12 ; (115) 

and the comparison with Eq. (46) reads: 

a 2 2 A 12 = ^e 2 31 A 3 ; (116) 

(,33 

in the limit of axisymmetric configurations, a 2 — > ai, Ai 2 — > An, the left- 
hand side of Eq. (116) takes the expression (e.g., Caimmi 1995): 



ta a ^ 12 = iU-4^f|i] ; (117) 
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and the combination of Eqs. (116), (117), yields: 



4^3 



5-44 + 4(l-4)Cn-Cs3 



3-2ei 



31 



3-2e 2 



31 



C: 



33 



e 2 i 



(118) 



which is an explicit expression of Wiegandt (1982a,b) criterion for bifurca- 
tion, with regard to homeoidally striated Jacobi ellipsoids in rigid rotation. 
In the limit of isotropical residual velocity distribution, (n = (22 = (33, 
Eq. (118) reduces to (Caimmi 1996a): 



£31^3 
A 1 



3 - 2e|i 
5-461! 



£21 = 1 



(119) 



as expected. 

The condition, defined by Eq. (118), has to be compared with its counter- 
part deduced in the current paper, under the same assumption i.e. C11 — C22 
also for triaxial configurations. The latter may be deduced from Eq. (77), 
using the relation (e.g., Caimmi 1966a): 



,. 2 A 2 -A x 

lim €r 



£21— >1 



-21 



1-4 




=31" 



A3-A1 
1-4, 



(120) 



and the result is again Eq. (118). Then Wiegandt (1982b) criterion for bifur- 
cation can be deduced from the current theory, provided isotropic residual 
velocity distribution is assumed along the equatorial plane, also for triaxial 
configurations. On the other hand, the above assumption is in contradiction 
with Eqs. (59) and then it cannot be accepted. A criterion for bifurcation, 
consistent with Eqs. (59), is expressed by Eq. (119). 

The origin of the discrepancy could be the following. Wiegandt (1982a) 
assumption (ii) implies Eq. (12) therein which, after integration, yields isotropic 
residual velocity distribution along the equatorial plane i.e. C11 — C22- On 
the other hand, an isotropic residual velocity distribution along the equato- 
rial plane is related to a zero-th order approximation in some special series 
developments [Marochnik 1967; Eq. (19) therein]. The additional require- 
ment of axial symmetry makes the above mentioned relation be exact, but 
in presence of triplanar symmetry it has to be considered as a zero-th order 
approximation. 
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In the light of the results discussed in the current subsection, the anal- 
ogy between the behaviour of collisional and collisionless self-gravitating flu- 
ids, subjected to the restrictions (Wiegandt 1982a,b): (i) rigid rotation; (ii) 
constant residual velocity on the boundary; (iii) absence of internal energy 
transports; appears to be complete. 
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